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1.0  INTRODUCTION 


Scientists  are  aware  and  the  general  public  is  becoming  more  and  more 
aware  that  life  on  earth  is  possible  only  because  the  atmospheric  ozone  layer 
exists.^  Ultraviolet  light  from  the  sun  is  absorbed  by  the  ozone  layer,  prevent- 
ing ultraviolet  rays  in  quantities  dangerous  to  human  existence  from  reaching  the 
2 

earth.  Were  the  ozone  layer  to  be  destroyed,  there  would  be  no  such  safety 
shield  protecting  life  on  earth. 

Medical  experts  may  have  linked  excessive  ultraviolet  light  exposure  to 

cancer  and  are  justifiably  concerned  that  even  partial  deterioration  of  the  ozone 

layer  would  allow  more  ultraviolet  light  to  penetrate  to  the  earth's  surface,  in- 

3 

creasing  in  number  the  victims  of  cancer. 

4 

We  know  that  ozone  is  very  unstable  and  that  it  disassociates  readily. 

The  ozone  "layer"  is  really  an  ozone  equilibrium,  a balance  between  ozone 

being  created  and  destroyed.  Atmospheric  contaminants  such  as  may  be  released 

from  aerosol  cans  or  jet  engine  exhaust  may  rise  high  enough  to  reach  the  ozone 

5 

layer.  If  they  do,  the  ozone  equilibrium  could  be  upset,  with  only  negative 
consequences. 

The  only  way  to  "keep  an  eye"  on  the  ozone  layer  is  to  measure  it,  and 
the  best  way  to  measure  it  at  altitudes  higher  than  30  km  is  via  sounding  rockets. 


An  ideal  way  to  measure  the  ozone  layer  is  to  incorporate  an  ozonesonde 


Jk 


SPACE  DATA  CORPORATION 

TEtFE , ARIZONA 


daily  from  facilities  around  the  world.  Such  a system,  if  feasible,  would  pro- 
vide a necessary  safety  measure  at  minimal  added  costs. 

The  detection  method  is  to  use  a sensor  disk  coated  with  a Rhodamine  B 
solution,  which,  when  contacted  by  ozone,  results  in  a light-emitting  chemical 
reaction  (a  chemiluminescent  reaction).  The  measurement  technique  is  to 
measure  the  intensity  of  the  light  emitted  due  to  the  ozone-Rhodamine  B reaction. 
Theoretically,  the  greater  the  intensity,  the  more  ozone  there  is. 

The  detection  system  must  be  designed  so  that  the  light  meter  measures 
only  light  emitted  due  to  the  chemiluminescence.  Extraneous  sunlight,  if 
visible  by  the  light  meter,  would  give  false  ozone  readings.  Therefore,  the 
ozonesonde  design  must  allow  flow  of  a representative  air  sample  while  simul- 
taneously blocking  out  stray  light.  One  proposed  design  is  shown  in  Figures  1. 1 
and  1.2. 

The  purpose  of  the  research  effort  described  in  this  report  was  to  define 
the  quantitative  relationship  between  the  intensity  of  the  emitted  light  and  the 
ozone  concentration  and  to  reveal  and  (if  possible)  solve  other  problems  intro- 
duced as  a result  of  the  research  effort. 

Consistency  and  repeatability  of  results  were  investigated  during  Sen- 
sitivity Tests.  Several  questions  arose  in  this  area,  none  of  which  precluded 
the  chemiluminescent  technique  for  measuring  ozone. 

Detector  output  versus  ozone  concentration  as  functions  of  flow  rate  and 
pressure  were  investigated,  though  on  a limited  basis  due  to  economic  constraints 


FIGURE  1 . 1 


EXHAUST 


FIGURE  1.2 

OZONESONDE  AIR  SAMPLING 
PROBE  CONFIGURATION 
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of  the  contract.  These  tests  revealed  the  importance  of  flow  rate  on  the  data 
reduction  technique.  Because  in-flight  flow  metering  is  too  expensive  to  be 
practical,  flow  rates  through  the  ozonesonde  ducts  must  be  calculated.  (See 
SDC  TM-1414.)  To  verify  flow  rate  calculations,  blowdown  tests  were  conducted. 

Inlet  losses  were  investigated,  but  they  were  found  to  be  insignificant 
within  the  range  of  tests  conducted. 

Finally,  vibration  tests  were  conducted  in  order  to  determine  whether 
vibrations  experienced  during  flight  would  be  sufficiently  forceful  to  cause  the 
sensor  material  to  flake  off,  resulting  in  decreased  sensitivity.  They  were  not. 

As  a result  of  these  tests,  the  conclusion  drawn  is  that  the  rocketborno 
ozonesonde  employing  the  chemiluminescent  technique  for  measuring  atmospheric 
ozone  is  feasible.  The  following  sections  describe  test  results  and  define  the 
atmospheric  ozone  data  reduction  technique. 


i 
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2.0  SENSITIVITY  TESTS 

2. 1 Introduction 

The  proposed  technique  for  measuring  atmospheric  ozone  incor- 
porates a photomultiplier  to  measure  the  light  emitted  as  a result  of  a chemi- 
luminescent reaction  between  ozone  and  Rhodamine  B,  a solution  of  which  is 
coated  on  the  ozonesonde  sensor  disk  surface.  Theoretically,  there  is  a meas- 
urable relationship  between  the  amount  of  ozone  reacting  with  Rhodamine  B 
and  the  intensity  of  the  resulting  emitted  light. 

The  purpose  of  the  Sensitivity  Tests  was  to  study  the  relationships  be- 
tween ozone  concentration  and  detector  (photomultiplier)  output.  Clearly, 
consistent  and  repeatable  results  (per  individual  disks)  are  mandatory  if  ozone 
measurements  are  to  be  accurate  and  valid.  To  facilitate  repeatability,  an  ozone 
generator  that  would  enable  known  amounts  of  ozone  to  be  created  at  repeatable 
rates  was  constructed.  Ozone  is  created  by  bombarding  a mixture  of  oxygen  and 
nitrogen  gases  with  ultraviolet  light  - as  is  done  in  the  atmosphere,  where  the 
sun  provides  the  ultraviolet  light  (UV)  source.  The  more  UV  there  is,  the  more 
ozone  is  generated,  and  the  relationship  is  a direct  proportion.  What  the  ozone 
generator  does  is  uncover/cover  a UV  lamp  at  a constant  rate,  such  that  more/less 
ozone  is  created.  A Dasibi  ozone  monitor  measures  the  amount  of  ozone  created, 
and  a strip  chart  recorder  can  be  (and  was)  used  to  record  the  Dasibi  ozone  con- 
centration measurements  and  the  photomultiplier  output.  Figure  2.1  is  a schematic 
of  the  test  setup. 
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FIG  LIRE  2.1,  TEST  SET  UP 
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2,2  Results  and  Observations 

Initial  tests  monitored  ozone  concentration  versus  detector  output. 
Results  of  these  tests  showed  that  consistent  and  repeatable  results  could  be  achieved, 
but  to  do  so  the  disks  had  to  be  sensitized.  Uhsensitized  disks  showed  consistent 
trends,  but  not  consistent  quantitative  results.  Sensitized  disks  showed  consistent 
trends  and  consistent,  repeatable,  and  predictable  quantitative  results.  To  say  disks 
are  "seraitized"  means  that  they  are  "warmed  up"  for  the  task  of  measuring  ozone. 
That  is,  ozone  concentration  versus  detector  output  stabilization  has  been  achieved 
in  sensitized  disks.  This  is  accomplished  merely  by  flushing  the  system  with 
ozonized  air  for  a specified  time. 

Further  testing  investigated  the  sensitization  phenomenon  to  answer  the 
questions; 

1)  If  exposure  to  ozone  sensitizes  disks,  does  zero  ozone  exposure 
desensitize  them? 

2)  If  so,  what  time  periods  are  involved  in  desensitization? 

3)  What  is  the  quantitative  relationship  between  a period  of  zero 
ozone  exposure  and  the  subsequent  and  immediate  ozone  flux 
required  to  resensitize  the  disks? 

4)  Specifically,  is  the  relationship  discovered  in  answering 
question  ^3  compatible  with  actual  flight  events?  That  is, 
during  a flight  there  will  be  a brief  period  of  zero  ozone 
exposure.  If  that  period  is  such  that  significant  disk 
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dosensitization  occurs,  how  much  subsequent  ozone  exposure 
is  required  to  resensitize  the  disk?  Is  the  amount  of  ozone 
present  at  high  altitudes  sufficient  in  quantity  to  accommodate 
disk  resensitization  at  flaw  rates  expected  during  flight? 

To  answer  these  questions,  the  same  test  setup  as  shown  in  Figure  2.1 
was  used.  The  test  procedure  was  to  sensitize  the  disks  by  flushing  with  ozone 
of  1 part  per  million  by  volume  (ppmV)  concentration  for  ten  minutes  at  1 liter/ 
minute  flow  rate.  A cycle  run  was  made  to  measure  ozone  concentration  versus 
photomultiplier  output  for  a sensitized  disk.  (During  initial  tests  of  this  kind 
two  more  cycles  were  conducted  and  the  subsequent,  repeatable  results  proved 
the  disks  to  be  sensitized.) 

The  disk  then  sat  idle  for  a specified  time,  after  which  cycle  runs  were 
made.  Data  taken  after  an  idle  period  were  compared  with  pre-idle  data,  and 
the  difference  showed  the  effect  of  zero  ozone  exposure  for  that  idle  period. 

Results  were  that  periods  of  zero  ozone  exposure  did  cause  disk  desen- 
sitization,  resulting  in  on  inconsistent  overshoot  phenomenon.  That  is,  desen- 
sitized disks  exhibited  a greater  photomultiplier  response  than  sensitized  disks. 
When  the  ozone  concentration  was  held  constant,  the  desensitized  disk's  photo- 
multiplier response  decreased  until  stabilization  (resensitization)  occurred.  The 
time  required  for  stabi lization  varied,  depending  upon  the  degree  of  desensitiza- 
tion and  the  ozone  flux  through  the  sensor  duct  during  resensititzation. 


-9- 


Figure  2.2  exhibits  the  overshoot  phenomenon  observed  in  desensitized 
disks.  The  particular  disk  shown  in  Figure  2.2  was  left  idle  for  fifteen  minutes. 
Figure  2.3  shows  the  same  disk's  photomultiplier  response  after  an  idle  period  of 
two  minutes. 

Comparison  of  Figures  2.2  and  2.3  shows  that  the  longer  the  idle  period 
of  zero  ozone  exposure,  the  stronger  the  disk  desensitization,  the  greater  the 
overshoot  and  the  longer  the  period  required  to  resensitize  the  disk. 

The  effect  of  overshoot  is  to  flatten  out  the  photomultiplier  response 
curve.  In  isolated  cases,  the  overshoot  effect  even  reversed  the  slope  of  the 
photomultiplier  response  curve!  Clearly,  if  photomultiplier  response  is  to  be 
used  to  show  the  amount  of  ozone  present,  the  overshoot  phenomenon  cannot 
be  present.  That  is,  disks  must  be  sensitized. 

Figure  2.4  shows  the  ozone  concentration  versus  photomultiplier  output 
for  a sensitized  disk.  Note  how,  after  an  initial  bend,  the  response  is  relatively 
linear  and  is  therefore  very  useful  for  measuring  ozone. 

Figure  2.5  shows  the  same  disk  after  an  idle  period  of  15  minutes.  Note 
how  the  photomultiplier  response  overshoots  the  reading  of  the  same  disk  when 
sensitized  for  the  same  early  ozone  concentrations.  As  the  disk  becomes  more 
exposed  to  ozone,  it  becomes  resensitized,  and  if  the  ozone  concentration  were 
held  constant,  the  photomultiplier  response  would  decrease  until  stable.  But  the 
ozone  concentration  is  increasing,  so  the  photomultiplier  response  is  also  in- 
creasing, although  its  rate  of  increase  is  decreasing  and  will  continue  to  decrease 
until  the  disk  is  resensitized. 
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Figure  2.6  shows  an  example  of  o drastic  overshoot  in  photomultiplier 
response.  Clearly,  this  response  is  unacceptable  for  ozone  measurements.  For 
this  case,  there  is  no  way  to  relate  photomultiplier  response  to  a specific  ozone 
concentration. 

When  the  disk  of  Figure  2.6  is  sensitized,  its  response  curve  is  stable, 
as  shown  in  Figure  2.7,  Clearly,  the  disk  - when  sensitized  - is  useful  for 
measuring  ozone. 

Thus, the  answers  to  the  first  two  questions  are; 

1)  Periods  of  zero  ozone  exposure  do  result  in  disk  desensitization, 
and 

2)  Idle  periods  as  brief  as  two  minutes  can  cause  measurable  disk 
desensitization. 

Since  the  time  between  liftoff  and  the  start  of  ozone  measurements 
is  on  the  order  of  two  minutes,  measurable  disk  desensitization  can  occur,  even 
if  disks  are  continuously  flushed  with  ozone  until  liftoff. 

Test  procedures  to  measure  resensitization  versus  desensitization  times 
were  as  follows; 

1)  Disks  were  fully  sensitized, 

2)  Disks  underwent  a specified  idle  period. 

3)  Disks  were  flushed  with  ozone  at  constant  concentrations, 
the  stabilization  was  recorded  and  the  stabilization  time 
was  measured. 

Results  of  these  tests  are  shown  in  Table  2. 1 . 


-15- 


FIGURE  2.6 


DESENSITIZED  DISK 


ALT  1 N.G  .d!r A5XLC  Q_\6tRSilO Q-t  EiLEKO-MXN ON 


Concentration 


. 

III. 

• • 

. ... . 

r 


i 


SPACE  DATA  CORPORATION 

TEMPE , ARIZONA 


More  meaningful  data  will  be  achieved  if  we  relate  resensitization  by 
ozone  flux  through  the  sensor  duct. 

From  Table  2.2,  a zero  ozone  exposure  period  of  two  minutes  requires 
an  ozone  flux  on  the  order  of  10^  molecules/second  for  two  minutes  in  order 
to  resensitize  the  disk. 

The  questions  ore:  what  is  the  ozone  flux  at  high  altitudes,  and  how 
long  will  it  take  to  resensitize  the  disk  with  such  an  ozone  flux? 

1 2 

The  calculated  ozone  flux  expected  at  altitude  is  on  the  order  of  10 
molecules/second.  Such  a small  ozone  flux  would  require  approximately  two 
hours  to  fully  resensitize  the  disk.  Such  a long  resensitization  time  seemingly 
should  preclude  testing  but  it  does  not  for  the  following  reasons: 

As  ozone  flux  decreases,  the  magnitude  of  the  overshoot  decreases. 

At  the  low  limit,  despite  the  overshoot,  the  photomultiplier  output  is  very 
nearly  equal  to  the  output  of  a stabilized  disk  for  the  same  concentration,  i.e., 
the  overshoot  is  negligible,  and  it  will  remain  negligible  until  greater  ozone 
number  density  measurements  need  to  be  made;  however,  when  the  ozone  con- 
centration is  high  enough  to  cause  the  overshoot  to  be  non-negligible,  the  re- 
sensitization times  for  those  higher  concentrations  are  shorter,  on  the  order  of 
two  minutes.  Since  the  descent  phase  of  the  flight  is  on  the  order  of  30  - 40 
minutes,  a two  minute  resensitization  time  is  acceptable. 

Also,  the  ozone  flux  does  not  remain  at  the  low  level  for  very  long.  It 
continually  increases  with  the  increasing  ozone  levels  present  at  lower  altitudes. 


-19- 


Observations 

The  sensitivity  tests  led  to  the  following  observations: 

1)  Exposure  to  ozone  sensitizes  disks. 

2)  Zero  ozone  exposure  desensitizes  disks. 

3)  Even  zero  ozone  exposure  periods  as  brief  as  two  minutes  result 
in  measurable  disk  desensitization. 

4)  The  effect  of  disk  desensitization  is  an  overshoot  in  photo- 
multiplier output,  i.e.,  there  is  a greater  photomultiplier 
response  for  a desensitized  disk  than  for  a sensitized  one.  The 
overshoot  decreases  with  decreasing  concentration. 

5)  Disks  desensitized  for  two  minutes  require  an  ozone  flux  of 
1014  molecules/second  for  two  minutes  in  order  to  resensitize 

the  disks. 

6)  During  flight  there  will  be  a period  of  zero  ozone  exposure 
of  approximately  two  minutes  (the  time  between  liftoff  and 
sonde  exposure  to  atmosphere);  therefore,  to  resensitize  disks 
in  flight,  an  ozone  flux  of  10^4  molecules/second  for  two 


minutes  is  required. 

There  is  not  enough  ozone  above  70  km  altitude  to  accommo- 
date the  ozone  flux  required  to  resensitize  flight  disks;  but, 
because  the  overshoot  in  photomultiplier  output  caused  by  de- 
sensitized disks  decreases  with  decreasing  ozone  concentration. 
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the  overshoot  is  low  for  low  ozone  concentrations.  Above  70 
km,  the  ozone  concentration  is  so  low  that  the  effect  of  over- 
shoot is  negligible. 

Conclusions 

1)  Disk  sensitivity  is  compatible  with  rocketborne  ozonesondes 
utilizing  the  chemiluminescent  technique  for  measuring  atmos- 
pheric ozone. 

2)  Disks  must  be  sensitized  prior  to  flight  by  flushing  with  1.0  ppm 
ozonized  air,  preferably  until  liftoff. 

3)  Disks  should  be  calibrated  as  close  to  liftoff  as  practical. 
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3.0  DETECTOR  OUTPUT  VERSUS  CONCENTRATION 

3. 1 Flow  Rate  Dependency 

Flow  rates  through  the  ozonesonde  sampling  duct  during  an  actual 
flight  are  expected  to  vary  between  6 and  300  liters  per  minute.  Flow  rates  of  6 
liters  per  minute  can  be  readily  achieved  in  the  laboratory.  Flow  rates  on  the 
order  of  300  liters  per  minute  require  expensive  vacuum  chamber/pump  systems 
that  were  not  available  within  the  financial  constraints  of  this  contract;  therefore, 
tests  measuring  detector  output  versus  concentration  as  a function  of  the  full  range 
of  flow  rates  expected  during  flight  were  not  conducted.  Flow  rate  dependency  of 
limited  scope  is  included  in  Section  3.2. 

3.2  Pressure  Dependency 

Low  pressure  tests  were  conducted  using  the  test  setup  schematicized 
in  Figure  3.1.  Essentially  what  was  done  was  to  sensitize  the  disk  by  flushing  with 
ozonized  air  for  approximately  twenty  minutes,  then  calibrating  at  ambient  pressure, 
evacuating  the  belljar  to  a low  pressure,  resensitizing  the  disk,  and  finally  cali- 
brating at  the  low  pressure. 

Figure  3.2  compares  results  at  ambient  pressure  with  those  at  low  pressure. 

While  the  system  was  set  up  at  low  pressure,  tests  were  conducted  to  show 
the  effect  of  small  changes  in  flow  rate  on  detector  output  at  low  pressure.  Those 
results  are  also  depicted  in  Figure  3.2. 

As  shown  by  Figure  3.2,  small  changes  in  flow  rate  result  in  measurable 
detector  output  differences.  During  flight,  large  flow  rate  differences  (on  the 
order  of  200  liters  per  minute)  are  expected. 
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FIGURE  3.1,  TEST  SETUP 


1 

* 

The  question  is,  what  effect  will  large  variations  in  flow  rate  have  on 
detector  output? 

To  answer  this  question  experimentally  requires  sophisticated  vacuum 
chamber  testing  beyond  the  scope  of  this  contract;  therefore,  we  must  answer 
the  question  analytically.  Appendix  A shows  how  the  detector  output  is  affected 
by  flaw  rate  and  altitude  effects. 
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4.0  BLOWDOWN  TESTS 

Previous  sections  have  shown  that  atmospheric  ozone  data  reduction 
depends  on  flow  rate  through  the  ozonesonde  duct.  Because  in-flight  flow- 
metering is  too  expensive,  flow  rate  versus  altitude  is  determined  by  calcula- 
tion. The  difference  between  the  theoretical,  calculated  flow  rote  and  that 
observed  experimentally  is  due  to  friction,  assuming  the  calculation  to  be 
otherwise  correct.  As  a result  of  the  blowdown  test,  the  calculated  friction 
factor  can  be  adjusted  to  make  the  calculated  flow  rate  agree  with  that 
measured  in  the  lab.  Friction  factors  are  independent  of  altitude;  therefore, 
if  the  flow  rate  calculation  technique  is  correct,  then  if  the  calculated  flow 
rate  agrees  with  the  observed  flow  rate  at  any  one  altitude  (the  test  altitude, 
for  example),  then  the  calculated  flow  rate  will  agree  with  the  observed  flow 
rate  at  any  and  all  other  altitudes. 

Flow  rates  through  two  ozonesonde  sampling  heads  of  different  designs 
were  measured  by  doing  blowdown  tests.  Each  head  was  mounted  (in  turn)  into 
a vacuum  chamber  such  that  the  ozonesonde  duct  entrance  was  at  ambient  pressure, 
and  the  exit  was  at  vacuum  chamber  pressure.  After  plugging  the  ozonesonde  duct 
entrance  and  sealing  all  sources  of  possible  leaks,  the  vacuum  chamber  was  evacu- 
ated to  a pressure  of  about  60  mm  of  mercury.  At  time  zero  the  duct  entrance  plug 
was  removed;  therefore,  air  flowed  through  the  sampling  duct  into  the  vacuum 
chamber.  By  timing  the  measured  pressure  increase  within  the  chamber,  the  flow 
rate  through  the  duct  was  calculated,  as  follows: 


-% 
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V_  dP 
P dt 


where  Fv 


V 

P 


dP 

dt 


= volumetric  flow  rote 

= vacuum  chamber  volume 

= vacuum  chamber  pressure 

= slope  of  chamber  pressure  versus  time  curve 

This  expression  is  derived  in  Appendix  B. 

Ozonesonde  blowdown  test  results  are  plotted  in  Figures  4.1  and  4.2. 
Vacuum  chamber  pressure  versus  time  is  shown  in  Figure  4. 1 , while  Figure  4.2 
shows  ozonesonde  sampling  duct  ^1  flow  rate  versus  back  pressure  ratio.  Back 
pressure  ratio  is  the  ratio  of  duct  back  pressure  to  duct  entrance  pressure. 

Using  techniques  defined  in  Appendix  C,  the  flow  rate  was  calculated. 
Differences  between  the  calculated  and  measured  flow  rates  for  ozonesonde  de- 
sign ^1  are  shown  in  Figure  4.2.  To  accurize  flow  rate  calculations  so  that  the 
flow  rate  agreed  with  the  measured  rate,  loss  factors  were  adjusted.  For  ozone- 
sonde head  design  *1,  the  loss  factor  was  changed  from  6.15  to  3.90.  For  ozone- 
sonde design  ^2,  the  loss  factor  was  estimated  to  be  3.00. 
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5.0  INLET  LOSSES  DETERMINATION 

Inlet  losses  were  determined  by  first  calibrating  using  the  setup 
schematicized  in  Figure  5.1,  then  calibrating  per  the  setup  shown  in  Figure 
5.2.  Comparison  of  results  yielded  the  inlet  losses. 

Within  the  range  of  tests  conducted  (flow  rates  varied  from  1 to  12 
liters  per  minute,  pressures  varied  from  10  to  100  millibars)  no  significant 
losses  were  measured. 


<5 AS  EXIT 


FIGURE  5.  1 , TEST  SETUP 
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6.0  VIBRATIONS  TEST 

A sensor  disk  was  mounted  on  an  ultrasonic  vibrator  and  shaken  for 
three  minutes  to  simulate  flight  conditions.  There  was  no  physical  evidence 
of  flaking  of  sensor  material  and  calibrations  run  prior  to  and  after  vibration 
testing  showed  no  decrease  in  sensitivity. 

Vibrations  expected  during  flight  will  not  cause  flaking  of  sensor 
material;  therefore,  vibrations  expected  during  flight  will  not  alter  disk 
sensitivity. 


4 
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7.0  THE  ATMOSPHERIC  OZONE  DATA  REDUCTION  TECHNIQUE 

Calibration  of  the  detection  system  will  be  in  the  form  of  ozone  con- 
centration (ppmV)  versus  detector  output.  The  calibration  curve  will  be  con- 
verted via  computer  to  ozone  number  density  versus  detector  output,  as  follows; 


<N°'eo, 


^ P o Peal  £ (10f6, 


where 


(No>ca. 


Because 


calibration  ozone  number  density,  molecules/cm' 

ozone  concentration,  ppmV 

Avogadro's  Number,  molecules/mole 

ozone  density  at  reference  pressure,  grams/liter 

ozone  molecular  weight,  grams/mole 

ozone  density  references  pressure,  inches  Hg 

calibration  ambient  pressure,  inches  Hg 

2 

conversion  factor,  liters/cm 
one  part  per  million 


6.02 x 1023, 

2.144, 

29.92, 


48,  and 


the  above  expression  reduces  to; 


(N  ) 

° cal 


8.9371  (10  ) C P^j  molecules/cm3 
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Thus,  the  calibration  curve  will  have  been  converted  to  ozone  number 
density  versus  detector  output. 

Flight  data  acquires  detector  output  versus  altitude,  which  derives  from 
the  calibration  curve  the  uncorrected  ozone  number  density  versus  altitude. 

Appendix  A shows  how  the  number  density  is  corrected  for  flow  rate 
and  altitude  effects  by  the  factor  ^ , where  Fv  = volumetric  flow  rate, 
and  ^ = mean  free  path.  Thus,  the  final  data  reduction  expression  is 


<N0>alt  = 

(No)*,,  x r where 

/Fv\ 

U;  cal 

(No>ol,  " 

ozone  number  density  at  altitude 

<N°>cal  " 

calibration  number  density  producing  the  same 
detector  output  observed  at  altitude 

■ 

volumetric  flow  rate  to  mean  free  path  ratio 
at  altitude 

o'  ^ 

Q_ 

II 

ratio  during  calibration 

By  reducing  the 

data  in  this  manner  for  a number  of  altitudes,  an 

atmospheric  ozone  profile  may  be  generated  and  plotted.  The  ozone  profile 
shown  in  Figure  7. 1 was  generated  using  the  data  reduction  technique  defined 
in  this  section.  The  computer  output  for  the  data  reduction  is  presented  in 
Appendix  D. 
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FIGURE  7.1,  ATMOSPHERIC  OZONF  PROFI  LE 


SPACE  DATA  CORPORATION 

PMOIMItt,  AtllOWA 


8.0  CONCLUSIONS 

Rocketborne  ozonasondes  can  use  the  chemiluminescent  technique  to 

measure  atmospheric  ozone.  To  do  so  most  efficiently  requires  sensitizing  the 

sensor  disk  by  continuously  flushing  with  ozonized  air  for  a short  period  prior 

to  and  up  to  liftoff.  Calibration  of  the  detection  system  is  required  as  close 

to  launch  as  is  practical.  Calibration  data  will  be  acquired  in  the  form  of 

ozone  concentration  (ppmV)  versus  detector  output.  Flight  data  will  yield 

detector  output,  from  which  may  be  derived  the  uncorrected  ozone  flux.  Final 

3 

data  reduction  via  computer  will  yield  ozone  concentration  (molecules  per  cm  ) 
versus  altitude.  A sample  computer  output  illustrating  reduced  ozone  data  is 
presented  in  Appendix  D. 
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APPENDIX  A 
THE  DATA  REDUCTION  EQUATION 


If  an  ozone  molecule  is  a distance  y from  the  sensor  plate,  then  the 
time  it  takes  the  molecule  to  diffuse  to  the  plate  is: 


L. 

ty 


Vy: 


* = Vy 


where 


then 


diffusion  velocity 


If  tx  is  the  time  in  which  the  elemental  volume  moves  down  the  tube, 
tx  = ~~  , where  Vx  = -j- 


where  Fy 
A 
(A 

Then: 

t 


x 

flow  rate 

cross-sectional  area  of  duct 

hW) 


L = _LA_ 
Vx  FV 


= critical  time.  This  is 


the  time  that  the  elemental  volume  is  in  contact  with  the  sensor  plate. 

Now  then,  if  ty  > tc,  then  the  ozone  molecule  a distance  y from  the 

A-l 


i 
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sensor  plate  will  not  be  sensed  (because  it  will  take  too  long  to  reach  the  plate). 
Said  another  way,  if  ty  < tc,  then  all  the  measurable  ozone  in  the  volume 
y Wdx  will  be  measured.  The  fraction  of  ozone  molecules  that  moves  toward 
sensor  plate  sensed  to  those  present  is  f,.,  where 


volume  of  ozone  measured  _ 6 y W dx  y 
volume  present  h W dx  6FT 


(There  is  a 1:6  chance  that  direction  of  Vy  will  be  toward  sensor.) 


Since  y 


such  that  t„<  t 
y — c 


Z_  ii. 


Thus 


LA  . LAVy 

FV  =* r * FV 

Vy 

= hW,  y < LhW  p^~ 

y _ LhW  Vy  = LW  Vv 

6K  6h  Fw  6 Fy 


i.e.,  the  fraction  of  ozone  measured  to  o zone  present  is; 
r . _ LW  V 


where  fx 


the  fraction 


sensor  length 


sensor  width 


diffusion  velocity 
flow  rate 


What  does  all  this  mean? 


A- 2 


When  we  measure  ozone,  we  get  a reading  0 which  indicates  how  much 
ozone  is  measured.  This  is  not  necessarily  the  amount  of  ozone  present.  In  fact. 


amount  of  ozone  present  = amount  of  ozone  measured 

hM 

i.e.,  Fm  = (ozone)  measured  =j>  (ozone)  present  = (ozone)  measured 

(ozone)  present 

Since  Fm  = l W Vy 

6 Fy 

(ozone  present)  = 6 [ (ozone)  measured]  Fy 

LW  Vy 


This  is,  essentially,  the  data  reduction  equation. 

What  we  must  do  now  is  relate  a calibration  to  altitude  measurements. 
During  calibration  we  hove: 

Fy 

(Oz)  present  = 6 x (Oz)  measured  x y 


vy  = k a , 

where  f{  = mean  free  path 

k = constant. 

This  says  that  the  diffusion  velocity  is  proportional  to  the  mean  free 
path;  the  greater  the  mean  free  path,  the  fewer  the  collisions  and  the  faster 
the  diffusion  rate. 

Thus,  during  calibration  we  have: 

lr^_\  _ 6 x (Oz)measured  x ^ x (®z)measured  * 

' ;pres®nt  LW  Vy  EwI7( 

and  we  can  plot:  0 versus  (Oz)  measured,  as  follows; 
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i 

< 

1 


Calibration 

0 

(Oz)m  ensured 

«1 

(No), 

02 

(n0)2 

®3 

• 

(nd)3 

• 

• 

en 

^o)n 

And  we  know;  Fy,  W,  and  /\  . 

During  flight  we  measure  0,  and  we  have: 


0 

(®z)mf»no  irftH 

eiF 

<No)1F 

02F 

(No>2F 

®3F 

• 

e 

(n°)3f 

• 

e 

• 

0nF 

(^o)nF 

Then;  f(°z)presentlfl;ght  6x  ((^measured!  flight 


(FV) 

LWl 


We  can  now  relate  flight  data  to  calibration,  as  follows: 
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t(°Z)  present]fMght  = LWk  A flight  ^ Unmeasured]  fl;ghf 

((°z) present] ca)  6 (FV>  «Oz)measured] 

LWk^col 

We  will  now  let  ( (Oi)measuracj]f|.ght  - , i.e.f  9f|jght  = 0caj 

L,  W,  and  k are  constant,  and  not  functions  of  altitude. 

Equation  reduces  to; 


f(°2)  present] f|;ght 

^flight  X cal 

'<0l)pr„.Mleo| 

/T  flight  (Fv)cq| 

or 

,(0Z)P™“">I  flight 

^°zVresent]ca|  x (^flight  Xcal 

(FV^cal  /T  flight 

This  is  the  data  reduction  equation. 

Thus,  the  data  reduction  equation: 

KOz) 

present  flight 

where 

I(0z)pr«®ntM,  x <FV>f|ight  x ^ ca| 

^FV^cal  ^ flight 

l(°z)present]ca| 

number  density  of  ozone  in  lab  that 
produced  the  same  0 reading  as  we  get 
at  altitude 

^vVlight 

= volumetric  flow  rate  at  altitude, 

calculated  by  computer  from  pressures, 
density,  fell  velocity,  etc. 

f\  flight 

= mean  free  path  at  altitude  (tabulated  in 

standard  atmosphere,  1976) 
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- calibration  flow  rate 

^ | = mean  free  path  at  calibration  altitude 


-% 
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APPENDIX  B 
BLOWDOWN  TEST  FLOW  RATE  DERIVATION 

Volumetric  flow  rate  will  be  determined  using  the  vacuum  chamber 
press  ire  versus  time  cirve,  the  known  (constant)  volume  of  the  vacuum  chamber, 
and  the  assumptions  that  the  air  acts  as  an  ideal  gas  and  the  temperature  remains 
constant. 

Thus,  the  eqmtion  of  state  for  an  ideal  gas  can  be  used: 


where 


e 

R 

By  definition,  p 


Thus, 

Since 


or. 


de 


dP, 


dP2 


dm 

dT 


where  F, 


Thus, 


v 

?FV 


fRT, 

Vacuum  chamber  pressure 

Density  of  air  in  chamber 
Universal  gas  constant 

tn  * where  m = mass  of  air  in  vacuum  chamber 

V V — chamber  volume 

J.  dm  (V  is  constant) 

V 

eRT, 

RT  d£  (R  is  constant,  T is  assumed  constant) 


RT  dm 
V~ 

V dP2  = 

RT  “ 


mass  flow  rote  = p F 

entering  chamber 


volumetric  flow  rate 

V dp2 

RT  d r 


B-l 


■ 


Sine* 
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p2 


V d P2 
■pT"  3T“ 

P RT  from  aquation  of  stot*. 


V dP2 

Pf  3 r 
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APPENDIX  C 

OZONESONDE  DUCT  FLOW  RATE  CALCULATION 


The  following  is  an  excerpt  from  SDC  TM-1414,  'Ozonesonde 
Duct  Flow  Rate  Calculation," 
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2.7  Finol  Flow  Rote  Equotion 

The  results  of  Sections  2. 1 through  2.7  are  conglomerated  to 
achieve  the  technique  for  calculating  the  flow  rate  through  the  ozonesonde 
sampling  duct  at  a given  altitude.  The  solution  is  as  follows: 

In  general. 


, ( 1 +2.507  V2; 

‘F  1 1 + 3 . 095  d /o' 


^lf  + lm  + 0 


•5  (*^r) 


pv\ 

T) 


Before  the  final  solution  can  be  achieved,  we  must  evaluate  the  compressibility 
)/  determine  whether  flow  is  laminar  or  turbulent  (and  thereby  establish  oC.  ), 
and  we  must  determine  whether  the  flow  is  in  the  continuum  or  transition  flow  regime. 

The  solution  thus  is  achieved  by  trial  and  error  and  is  best  achieved  via 
computer. 

The  quantities  A,  AeX;j.,  Lp,  L^,  d,  and  Kp  are  functions  of  the  duct  corv- 
figuration  and  are  constant  throughout  the  flight. 

The  quantities  Poo  , f ^ , and  /{  are  functions  of  the  altitude  in  question. 

Py  and  Pj  are  calculated  as  functions  of  Pqj  and  sonde  fall  velocity,  Uq, 
as  given  in  Section  2,1. 

©C  and  ^ are  functions  of  the  flow  rate.  Therefore,  their  values  must 
be  assumed  to  achieve  the  first  approximation.  The  first  approximation  will  give  us 
an  idea  of  whether  the  flow  is  laminar  or  turbulent  and  of  what  the  compressibility  may  be. 


I 
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The  solution  technique  is  outlined  as  follows; 


1) 

2) 


3) 

4) 

5) 


6) 


Using  the  methods  of  Section  2. 1 , calculate  p^“ 
Evaulate  Q , as  follows; 


If 

If 


»V  > .84,  $ =1. 
PT 


If 


PV  < .63,  B =0.7936. 

PT  r 

.63  < PV  < .84,  (ft  = 4.2404/ -.63\2+0. 

PT  ^ W ) 


7936 


LetO<_  = l. 

Determine  P®  and  for  the  altitude  in  question. 

Calculate 


= F 


Calculate  M^,  the  flow  mach  number  at  entrance 


I 


0.91  a^A. 


in 


where  F.  is  the  quantity  evaluated  in  Step  5), 

I 

= speed  of  sound  at  altitude  in  question. 


Joo 


A-  = duct  entrance  area 

in 


7)  If  M <:  1 and  if  ft  = 1,  ca  Iculate  Reynolds  number  per  Step  8). 


If  Mj  < 1 and  ^3  ^ 1,  then  changey^  to  1 and  recalculate  F.. 
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®)  If  Mj  > 1 (or  if  M j -<  ] and  = 1), 

R„  ~ Fj  fco  d = Reynolds  Number 

A J{co 

9)  If  Re  S:  800,  flow  is  turbulent  and  assumption  that 

©<-=  1 was  valid 

If  Re  ■<  800,  c4  = 800  , and  you  must  recalculate  F . 

R i 

e 

(i.e.,  return  to  Step  5)  and  redo  5)  through  8) 

• t 

10)  Re  = new  Reynolds  Number.  If  Re  is  within  10%  of  the  old 
Reynolds  number,  oCapproximation  is  close. 

N)  Calculate  the  ratio  d/2^  , the  ratio  used  to  determine  the 

flow  regime. 

If  ^/2({  = > 100,  continuum  flow  exists,  and  Fy  = Fj. 

If  d/2  R = — 0. 1,  purely  free  molecular  flow  exists, 

and  Fy  = Kp  ^ 1 - jV.  j . 

1^)  If  0. 1 c d c 100,  transitional  flow  exists,  and 

tt 

FV  = f|*ZKf  ,„h«r. 

Z = 1 + 2.507  d/2j 

T+T. 095d^X 


C-4 
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flow  rate  at  the  altitude  in  quesrion 
duct  cross-sectional  area  (in  sensing  chamber) 
duct  exit  area 
duct  entrance  area 

duct  hydraulic  radius  (in  sensing  chamber) 
friction  factor 
momentum  loss  factor 
duct  conductance 
ambient  pressure  at  altitude 
ambient  air  density  at  altitude 
viscosity  of  air  at  altitude 
mean  free  path  at  altitude 
speed  of  sound  at  altitude 
sonde  fall  velocity  at  altitude 
Venturi  pressure 

total  pressure  incident  to  duct  = P + — 

o>  2 

friction  parameter 
compressibility  parameter 
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Calculation  of  Py,  the  Venturi  (Vesture 


Assume  isentropic  flow  (no  losses).  Flow  characteristics  can  be  determined 
using  isentropic  flow  tables. 

• If  A*  >•  then  M - 1 in  throat. 


Thus,  for  all  ARAT1  such  that  ^ < ^4 


A>  A i 


, M = 1 at  throat. 


(A*  > A^^jf  * .<  4_ 

^ A Athroat  •) 


from  the  Appendix,  A throQt  = .2167  in2  and  A4  = .7854  in2.  Thus, 

Athroat  = A2  = 0.2167  in2,  A4  _ 0 7854  = 3 6244 

A2  0.2167 

For  all  ARAT1 , such  that  ~ -c  3.6244,  M = 1 at  throat. 

This  corresponds  to  —1  of  -<.982.  Thus,  for  all  ~ < .982/ 

fy  = 0.5283  Pp.  Since  P4  = P^,  , M = 1 at  throat  for  ^ ^ .982 

pT 

For  __4  >.982,  >3.6244,  subsonic  flow  exists  in  throat,  and  P 

Pt  A*  V 

is  found  using  isentropic  flow  tables. 

Summary  of  calculation  of  Py  using  isentropic  flow  tables* 
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SPACE  DATA  CORPORATION 

TDPC.  ARIZONA 


Calculate  PRATll  = If  PRATll  c .982, 

PT 


Pv  - .5283  PT. 

Look  up  PRATll  in  isentropic  flow  tables.  Find  corresponding 
ARAT11  (ARAT11  =A 4/A*). 

Calculate  ARAT22,  where  ARAT11  = ARAT1 1/3.6244. 

Look  up  ARAT22  in  isentropic  flow  tables  (subsonic  section). 
Find  corresponding  PRAT22. 


PRAT22  = ^ . So  Pv  = PRAT22*Pt. 
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